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Abstract 

From dumbbells to FCC crystals, we study the self-assembly pathway of amphiphatic, spherical 
colloidal particles as a function of the size of the hydrophobic region using molecular dynamics 
simulations. Specifically, we analyze how local inter-particle interactions correlate to the final self- 
assembled aggregate and how they affect the dynamical pathway of structure formation. We present 
a detailed diagram separating the many phases that we find for different sizes of the hydrophobic 
area, and uncover a narrow region where particles self-assemble into hollow, faceted cages that 
could potentially find interesting engineering applications. 
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INTRODUCTION 



Spontaneous assembly of components into large ordered aggregates is a ubiquitous phe- 
nomenon in nature, and is observed across all length scales. Aggregation of proteins into 
functional nanomachines [1], formation of viral capsids [21 [3], packing of phospholipids into 
biological membranes [1] and assembly of colloidal particles into macroscopic crystals [5] are 
just a few manifestations of this fundamental process. Because of advances in particle syn- 
thesis [U [TJ [HI ini [To] , it is now possible to produce colloidal particles that are anisotropic both 
in shape and surface chemistry, thus providing an unlimited variety of building blocks that 
can spontaneously assemble into an unprecedented number of structures holding promise for 
the development of materials with novel functional, mechanical, and optical properties. 

Although the generic features of particle aggregation can be described, at least phe- 
nomenologically, in terms of simple thermodynamic arguments [HI [121 (131 E] , the details 
of the process are far from being understood. In fact, for self-assembly to take place, a very 
delicate balance between entropic and energetic contributions, coupled to a precise geomet- 
ric character of the components, must be satisfied. In general, self-assembly is not to be 
expected unless a careful design of the building blocks has been performed beforehand, and 
this has inspired a large body of work dedicated to gain insight into how the geometry of the 
interparticle interactions and the shape of the particles themselves determine the dynami- 
cal pathway and structure they aggregate into. The ultimate goal is to be able to gather 
a sufficient understanding of the forward self-assembly process to then be able to develop 
tools that will allow us to tailor interparticle interactions to target desired structures. With 
this in mind, there is a clear necessity to explore a new dimension in the classical tempera- 
ture/concentration phase diagrams: the geometry of the interactions (see reference [15] for 
a recent perspective on the subject and a review of the relevant literature.) 

Recently, the self-assembly pathway of Janus particles (spherical particles with one hemi- 
sphere hydrophobic and the other hydrophilic) has been described in some detail [TB] both 
experimentally and numerically. Janus particles are important because they represent what 
is probably the simplest model where one can explicitly study and ask fundamental ques- 
tions about the role of surface anisotropy in colloidal aggregation. What was found is a 
rich behavior in terms of both the dynamics of self-assembly, and the final structure of the 
aggregates. Specifically, by tuning the repulsion between the hydrophilic hemispheres, one 
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can drive the system into two distinct phases. In one phase particles form a gas of small 
clusters composed of 4 to 13 particles, in the other phase, at large salt concentration, parti- 
cles organize into long, branched, worm-like structures formed by cooperative fusion of the 
small clusters. Such particles are a particular type of a wider class of patchy particles, the 
self-assembly of a wide variety of which has been studied previously p!7l [18] . 

Inspired by these experimental results, in this paper we go one step forward and explore 
dynamics and structure formation of anisotropic amphiphatic particles. While the dividing 
surface between hydrophobic and hydrophilic regions in Janus particles is set at the equator, 
we shall consider particles in which this boundary is located at an arbitrary latitude on the 
particle's surface, and we systematically study how, at a fixed concentration, these particles 
self-assemble depending on the size of the hydrophobic region. We thus explore a key 
"anisotropy dimension" proposed in [15] as related to the self-assembly of these particles. 

This system is quite interesting because one can easily predict the formation of a wealth 
of different structures as a function of the location of the dividing surface. For instance, we 
know that when the hydrophobic region covers a small area, particles can only assemble into 
dumbbells. We know that when the dividing surface is placed to slightly larger values, so that 
more particles can share the same hydrophobic surface, particles condense into small stable 
clusters. We also know the structures resulting from self-assembly of Janus particles, and of 
course, in the limit of full coverage, we recover an isotropic potential which is known to lead 
to the formation of an FCC crystal. Using this simple system we can study systematically 
the formation of aggregates whose structure ranges from zero dimensions (dumbbells and 
meso-particles) to three dimensions (FCC crystals), and we can analyze how the specificity 
of the local inter-particle interactions correlates to the final self-assembled structure and its 
dynamics. 

METHODS 

We model the amphiphatic character of the particles via an interaction potential that 
depends on both the separation between particle surfaces and the angle between particle 
axes, so that a precise shape and extent of the interaction may be defined. Our choice of 
the interaction potential is inspired by the model introduced in reference [16] that has been 
used to analyze actual experimental data of Janus Particles. 
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The potential presented in |T6] was specifically tailored to describe the physical prop- 
erties of Janus particles at different salt concentrations, ps- At low and intermediate salt 
concentrations, the repulsion between the charged hemispheres constrains particles to inter- 
act head-on, with angular deviations that depend in a nontrivial way on p^. However, the 
more interesting structure formation (worm-like clusters) appears at large salt concentra- 
tions, where the role of electrostatic interactions becomes negligible. In this limit, particles 
can freely rotate, once in contact, within the boundaries of the hydrophobic regions. These 
regions can therefore be appropriately described by a short-ranged, isotropic attractive po- 
tential that acts within the boundaries of that region. Similarly, the hydrophilic region is 
well-characterized by a simple short-ranged repulsive interaction. The potential used in our 
paper reflects this phenomenological behavior and has the following form: 



where r is the distance between particles, 6i is the angle between the axis of particle 1 and 
the axis between particle centers, and 62 is the analogous angle for particle 2. Vy.ep{r) is a 
symmetric repulsive interaction that accounts for the particle excluded volume, and has the 
form of a shifted-truncated Lennard- Jones potential: 



V^tt(^) accounts for the attraction between the hydrophobic surfaces on the particles. If 
Tg = \r — a\is the distance between the particle surfaces. 



and it extends up to r = 1.5cr. Finally, {0) is a smooth step function that modulates the 
angular dependence of the potential, and is equal to 1 within the region 6 < 6'jnax and decays 
to zero following the expression cos^(7r(^^ — O^s.^) / {26taii)) at the tail of the angular range, 
i.e when 6'max < ^ < ^max + Qtaii (witli Otaii = 10°). This particular value of Otaii has been 



illustration of {0). 

The parameters set in our model are compatible with colloidal particles of size (200 — 
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V{r,e,,e2) = Kep(r) + Ktt(r)0(^i)0(^2) , 





(1) 



selected to generate a sufficiently smooth potential at the Janus interface. See Fig. [T] for an 








FIG. 1: (Color online) Sketch of the angular dependence of the inter-particle potential (j){0). The 
dark side represents the hydrophobic region and the light side represents the hydrophilic. 

400)nm in aqueous solution kept at a salt concentration sufficiently large to screen the 
electrostatic repulsion between the hydrophilic regions of any two particles |16j . 

The system is evolved using Molecular Dynamics simulations with a Langevin thermostat 
at constant room temperature, T, in a cubic box with periodic boundary conditions. Our 
system contains = 10^ particles kept at a constant volume fraction (p = 0.01. We have 
chosen this concentration because it is comparable with those used in experimental studies 
on Janus Particles [16] , so that our work could have a grounding experimental reference that 
our results could be compared to for ^^ax — 90°. Each simulation runs for a minimum of 
10^ steps with a time step 6t = 0.001. All quantities in this paper are expressed in standard 
dimensionless units. 



RESULTS AND DISCUSSION 

Our goal is to understand how particles condense into stable three-dimensional aggregates 
via the process of self-assembly, and how the specificity of the geometry of the interaction is 
reflected in the final structure. Figure 2 reports one of the main results of our simulations. 
It shows a diagram indicating the self-assembly lines separating the structures obtained for 
different values of 6'max and e, with a typical resolution of one degree for 6'max, and O.I/cb^^ 
for the binding energy. 

As expected, a rich variety of structures arises depending on the position of the dividing 
surface ^max, and particles' binding energy, e. Notice that, consistent with recent numerical 
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simulations on self-assembly of viral capsids [19], and chaperonins [20], intermediate ordered 
aggregates are extremely sensitive to the size of the hydrophobic region and self-assembly 
occurs in a very narrow region. At low binding energy the system is in a gas state. For small 
values of 6'max and moderate values of e (Fig. |2](a)) particles aggregate into small clusters 
of 4 to 13 particles (meso-particles) including icosohedral structures like those seen in [2T] . 
For e 

max ~ 90°, self-assembly yields worm-like extended structures (Fig. |2|(b) ) , as has been 
observed previously for Janus spheres [16]. As the angular location of the hydrophobic 
region increases from 90° to 180° we find, in order of appearance, self-connected worm-like 
structures (Fig. [2]^c)), flat-crystalline bilayers (Fig. [2|^d)), faceted hollow cages (Fig. |2](e)) 
(similar to structures seen in the self-assembly of cone-shaped particles [221 EH], yet not 
restricted to specific magic numbers of particles.), amorphous fluid blobs (Fig. |2]^f)), and 
finally FCC crystals (Fig. |2|g)). 

Apart from the possible relevance of phases (d) and (e) for practical applications, we 
want to point out that none of these structures arises following a simple particle-by-particle 
growth mechanism, but rather by two-to-three step hierarchical self-assembly. The first step 
typically involves the formation of more- or less-structured meso-particles (depending on 
^max)- Next, small clusters organize into either extended worm-like aggregates, which then 
coalesce or deform to produce the final structure, or into larger fluid clusters that, once 
beyond some threshold size, spontaneously organize into structured aggregates. 

It is of particular interest to discuss in some detail the dynamical pathway of structure 
formation relative to phases (c), (d) and (e) as they all form via a complex three step 
mechanism (the dynamical pathway leading to structures in region (b) is identical to what 
was found in reference [16] and we refer the reader to that paper for a thorough description.) 
Surprisingly, the common precursor to all of them is the worm-like structure stable at ^max ~ 
90°. Self-connected worm-like aggregates are a consequence of the improved flexibility of 
the worm-like structures. As ^max increases, so does the ability of particles to rotate about 
their axes. The net result is that the branching ends of the clusters begin to connect, thus 
forming topologically nontrivial aggregates. 

To quantify the statistical difference between the aggregates found in regions (b) and 
(c), we measure their average radius of gyration Rq = (l/A'cX]j=i ^ (^"'''cm)^ >) as 
a function of cluster size. Results are shown in Fig. [3| and clearly indicate how as 6'max 
increases, the typical size of the aggregate decreases accordingly, resulting in more compact 
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structures. The low statistics for large clusters prevent us from making any meaningful 
estimate of the chains size exponent. We also measure the angular probability distribution 
function, P(cos(a)), between neighboring particles in clusters containing at least 50 parti- 
cles. Fig. |4] shows P(cos(a)) in regions (b), (c) and (d). P(cos(a)) shows a clear double 
peaked shape in region (d), where the aggregates assume a planar bilayer configuration, but 
the more intriguing difference is between region (b) and (c), highlighted in the insert of the 
figure. Region (b) is characterized by a peak at cos(a;) ~ — 1, and less distinguishable peaks 
at cos(a) ~ ±0.74 and cos(a) ~ —0.5; overall there is a large probability for all possible 
orientations. These data suggest that branched worm-like aggregates have a roughly circular 
cross-section with particles oriented in a disordered fashion, but for a slight preference for a 
few selected angles (remnants of the meso-particle structures they self-assembled from), and 
an anti-parallel neighbor opposite to most particles. In contrast, region (c) is characterized 
by a more distinct double-peaked function, with each peak close to ±1. This is consistent 
with a cross section that has flattened with respect to structures in region (b), and indicates 
that each branch of the aggregate is acquiring bilayer-like features. The dashed line separat- 
ing regions (b) and (c) shows where in the diagram the strings acquire sufficient flexibility 
to begin to form complex self-connected aggregates. 

This observation is quite interesting when related to the dynamic pathway leading to 
phase (d). In fact, bilayers are formed by either coalescence of co-planar loops, as shown 
in Fig. [5} or by branch alignment. This transition occurs at ^^^ax — HO", and it is driven 
by the large energy gain attained by the aggregates when each particle surrounds itself 
with the suddenly increased number of neighbors compatible with the enlarged hydrophobic 
region. This transition is quite sharp and represents a beautiful example of how a very 
small perturbation of the geometry of the local interaction can lead to completely different 
macroscopic structures. 

Finally, formation of finite-size faceted capsids does not occur by self-assembly of mis- 
oriented, disjoined bilayers, but via a mechanism which, once again, involves multiple steps. 
As we increase the size of the particles' hydrophobic region, short worm-like clusters, which 
are now very flexible, immediately fold onto themselves to form small, amorphous fluid 
blobs. We find that these blobs tend to remain fluid-like, whereas larger ones morph into 
faceted cages via a mechanism similar to that described for the formation of bilayers. Fusion 
of fluid blobs, as illustrated in Fig. [6} is the main mechanism through which large blobs. 
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which eventually turn into faceted cages, are generated. The dashed line between regions (d) 
and (e) indicates the onset of cage formation; however, a non-negligible number of planar 
aggregates are found to coexist with the faceted cages in region (e). It is worth stressing 
that these faceted cages can be considered as colloidal analogs of lipid vesicles; they are 
hollow, their inner walls are hydrophilic, and one could in principle consider using them 
as possible drug carriers, with edges and corners presenting convenient locations for outer 
surface tagging. 

A statistical analysis of our data correlating size and structure for non-planar aggregates 
in region (e) indicates a clear preference for large clusters to develop into faceted hollow 
cages (see Fig. [?]). What sets the onset cluster size for this transformation is a complex 
compromise between the geometric constraints imposed by the interparticle potential; the 
energy gain to close-pack particles in a bilayer, which grows with the number of particles A^; 
the energy cost for sides and corners, which have on average fewer neighbors than in a fluids 
state and have an energy cost which grows as N^^'^ and respectively; and finally the 
entropy loss due to particle ordering. Clearly, as increases, at sufficiently large binding 
energy, planar configurations become the most stable, and this results in surface faceting. 

When 6'max increases to even larger values, the geometric constraints imposed by the 
inter-particle potential become less restrictive, and planar bilayer configurations become 
less stable. Phase (f) is characterized by fiuid, amorphous blobs which remain fiuid at all 
sizes. Large clusters, formed by the smooth fusion of smaller ones, tend to contain smaller 
sub-clusters in their interior. This is the first sign indicating that aggregates begin to acquire 
a three-dimensional character, which eventually leads to the formation of clusters with FCC 
order, as inner and outer clusters begin to interact with each other. The location of the fiuid 
to FCC transition is at ^max — 135°, and was found by performing a structural analysis of the 
aggregates. Following [211 125], we identified particles whose local orientation is compatible 
to FCC ordering via a local bond order parameter based on spherical harmonics. Given a 
particle z, we consider 

Qnmii) = — — V F6m(ry) , (2) 
where j runs over the Nb{i) neighbors of particle i, from which a rotationally invariant order 
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parameter correlating the orientation of neighboring particles i and j can be defined as 

6 

q6(^) • qeO') = Yl Qemi^) ■ QlmU) (3) 

m=— 6 

Figure [s] shows the degree of crystallinity (O) as a function of 9^s.x across the angular range 
6'max £ [120 : 180]. (O) is obtained by first averaging qQ^i) ■ qeO) over the neighbors of each 
particle i, and then by taking the average over all crystalline particles in a cluster. Clearly, 
once crystalline particles are formed, their degree of order in a face-centered cubic crystal 
structure is independent of ^max- The explanation of this behavior is purely geometrical. 
In fact, for sufficiently large O^g^^, the hydrophilic area on each particle becomes so small 
that can be positioned among the 12 particles' contact points resulting from an FCC crystal 
without affecting its structure. Simple geometric considerations suggest that the onset value 
should occur for an angular span of the hydrophobic region larger than 150°. This value 
is compatible with our result 6'max — 135° because of the extra tail of 10° in our definition 
of the potential angular dependence. A careful analysis of the model for decreasing values 
of ^taih not shown here, does indeed result in a systematic shift of the crystallization onset 
6'max to larger values. 

At last, it is important to point out that, once again, the pathway leading to the for- 
mation of crystalline aggregates occurs in a two-step fashion. Particles first condense into 
a large, fiuid aggregate, and then crystallize from within the cluster via a standard nucle- 
ation process. Figure |9] shows in the same plot the size dependence of the whole aggregate 
over time, versus the number of crystalline particles in it. Clearly, crystal formation begins 
long after the aggregate is formed, and the large fiuctuations in the initial stages of crystal 
growth show the typical signature of crystal nucleation. These results are compatible with 
nucleation studies of isotropic colloidal particles interacting via a short-range attractive po- 
tential first observed in |26j, and highlight the crucial role played by meta-stable phases in 
the dynamics of crystal growth. In fact, we believe that this is the over-arching physics 
behind the rich dynamical phenomenology we find throughout this paper. Self-assembly 
proceeds, as predicted by Ostwald's step rule in the context of crystal nucleation [27], in a 
step-wise fashion that accounts for the complex free energy landscape containing multiple 
meta-stable states. Although we have not looked at the stability of the different phases 
found in our simulations, planar and vesicular structures are also being observed in the 
study of the equilibrium properties of a model system similar to ours [28]. It is also worth 
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mentioning that we expect the precise location of the phase boundaries to be somewhat 
sensitive to the particular choice of the angular potential. Unfortunately, this is very hard 
to characterize experimentally near the Janus interface and there is not a unique way of 
modeling that boundary. Our phase diagram is therefore intended to serve mostly as a 
guide for experimentalists. 

CONCLUSIONS 

In this paper we used molecular dynamics simulations to study the self-assembly pathway 
of spherical amphiphatic colloidal particles. We uncovered a wealth of different aggregates 
whose structures span the three dimensional spectrum. Specifically, depending on the size 
ratio between the hydrophilic and hydrophobic regions, particles self-assemble into small 
micellar clusters, worm-like structures, planar bilayers, faceted and fluid cages, and finally 
FCC crystals. We described the hierarchical self-assembly pathway leading to most of these 
structures and discussed their connection to the geometry of local inter-particle interactions. 
Finally, we made precise predictions for the formation of hollow amphiphatic cages. 

Although the morphology of some of our aggregates can be predicted by simple geometric 
considerations, the dynamics leading to their formation is far less trivial, and may play 
a crucial role in the efficiency of the self-assembly process. We believe that, apart from 
trivial cases, any procedure attempting to design inter-particle interactions to target specific 
structures could greatly benefit from taking into account the dynamics of structure formation 
in the design process. 
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FIG. 2: (Color online) Self-assembly diagram of amphiphatic colloidal particles as a function of 
binding energy e and size of hydrophobic region ^max* Region (a) is populated by small micellar- 
like clusters containing 4 to 13 particles. Region (b) contains branched worm-like aggregates. In 
region (c) we find self-connected worm-like aggregates. Region (d) delimits flat hexagonal bilayers. 
In region (e) we find faceted hollow cages. Region (f) is populated by fluid amorphous aggregates. 
Finally, in region (g) we find large clusters with FCC order. 
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FIG. 4: (Color online) Angular probability distribution function, P(cos(a)), between neighboring 
particles in aggregates containing at least 50 particles in region (b||0max = 91°)5 (c||^max = 95°), 
and (d||^max = 101°), of the phase diagram. The insert highlights P(cos(a)) in region (c) and (d). 
The data are averaged over 10 different clusters of each kind. 
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FIG. 5: (Color online) Sequence of three snap-shots from our shuulations showing the mechanism 
of bilayer formation via coalescence of three coplanar loops. 
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FIG. 6: (Color online) Sequence of six snap-shots (1^6) from our simulations showing faceted 
cage formation via fast folding of short worm- like clusters (1-2) and subsequent fusion of fluid 
blobs (3-6). 
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FIG. 7: (Color online) Probability distribution function P(n) as a function of cluster size n for 
fluid and faceted aggregates in region (e) of the phase diagram. 
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FIG. 8: Degree of crystallinity of self-assembled aggregates as a function of hydrophobic area 9i 
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FIG. 9: (Color online) Linear-Log plot of the typical aggregate size n as a function of time. The 
dotted line shows the total number of particles in the aggregate, while the continuous line shows 
how many of those particles are tagged as crystal-like. These data were collected at ^max = 160°. 
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